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SPHEROID AND DISK SIGNATURES IN GALAXY BULGES 
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Abstract. Recent progress on the structure and dynamics of bulges 
is reviewed. Those aspects that link galaxy bulges either to oblate 
C^") ' spheroids akin to elliptical galaxies or to rapidly-rotating, flattened 

systems more nearly resembling the products of disk internal transfor- 
mations, are highlighted. The analysis of surface brightness profiles 
derived from HST data is reviewed to show that unresolved nuclear 
components detected by HST have biased the determination of surface 
brightness profiles obtained from the ground; r 1 / 4 profiles are virtu- 
ally nonexistent in galaxy bulges. Predictions from accretion N-body 
models on the shape of surface brightness profiles are discussed. The 
' position of bulges on the Fundamental Plane (FP) of elliptical galaxies 

is examined to infer clues on bulge population ages and bulge dynamical 
structure. Kinematic diagnostics on the internal dynamics of bulges are 
examined. Finally, a new approach to the kinematic analysis of galax- 
ies, based on the use of synthetic spectra of single stellar population 
O ' models instead of the standard stellar templates, is presented. 
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1 Introduction 

In mature disk galaxies of early to intermediate Hubble type, the central regions 
bulge above the disk: in edge-on views, the isophotes defining the disk bend up 
d when approaching the galaxy center. Such central bulges have other distinct prop- 
erties; foremost, they have a smooth appearance in edge-on views, implying little 
dust and no current star formation activity outside of the plane of the disk. This 
motivates the notion of bulges being a separate galaxian component, with struc- 
ture, dynamics, populations and formation history distinct from those of the disk. 
The prominence of the central bulge is an important discriminant of the Hubble 
sequence, the other being the morphology of the spiral arms (Sandage & Bedke 
1994; see Gilmore 1999). 

The high luminosity densities of bulges make them accessible to observation 
even at high redshift. On the other hand, their small sizes relative to those of the 
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disk (ratio of scale sizes r e /h ~ 0.1; Courteau et al. 1996; Graham 2001, hereafter 
G01; MacArthur et al. 2003), make the study of the inner structure of bulges 
difficult; in small bulges, r e is comparable to the scale height of the disk. 

Work in the eighties led to the view that galaxy bulges are in most ways 
identical to elliptical galaxies that acquired a stellar disk around them. These 
similarities supported bulge formation models based on the monolithic collapse 
paradigm, or on mergers as explored by, e.g., Kauffmann (1996). During the past 
decade, observational results and theoretical modeling have provided new clues to 
re-address the basic questions on the nature of bulges: How similar are bulges to 
elliptical galaxies? How old are bulges? Is bulge formation prior to or after disk 
formation? Do bulges grow such that galaxies evolve along the Hubble sequence? 

Recent extensive reviews on bulges are found in Wyse et al. (1997) and Carollo 
et al. (1999). In this paper bulge properties are reviewed to highlight evidence for 
the spheroid vs. disky nature of bulges. Since the Kormendy (1993) review, this 
field has grown, mostly due to HST imaging work and to new kinematic studies. 

2 Definitions of bulges 

In §^we gave the original definition of a bulge: a smooth light distribution that 
protrudes from the central part of the disk in highly-inclined galaxies. This defini- 
tion emphasizes the old age of the bulge stars, the lack of both star formation and 
patchy dust. To implement this definition in practice, we look for a central region 
with lower ellipticity than that of the disk. This isophotal criterion was used by 
Kent (1986) and by Andredakis et al. (1995, hereafter APB95). 

Component analysis of galaxy structure (Kormendy 1977) introduced a pow- 
erful, easy-to-apply methodology to characterize bulges via fitting analytic surface 
brightness profile models such as the r 1 / 4 model. This led to a surface-brightness 
definition of bulges as the extra central light above the inward extrapolation of the 
constant-scale length exponential disk. The STScI workshop on bulge formation 
(Carollo, Ferguson & Wyse 1999) took this as the operational definition of a bulge. 
Widely used today, this definition might not be ideal, as it leads to assigning to 
the bulge any inner brightening of the disk. 

The co-existence of two definitions (isophote-based and surface brightness- 
based) for bulges is less than ideal. While the isophotal definition associates bulges 
with spheroidal structures, and suggests formation models similar to those of other 
galaxian spheroids such as elliptical galaxies, the surface-brightness criterion al- 
lows one to find bulges in late-type dwarf galaxies, as long as the surface brightness 
profile shows a central excess with respect to an exponential law. Most likely, the 
latter bulges do not share the same formation mechanism as the large bulges of 
early- type galaxies. 

3 Early results 

At the end of the eighties, bulges were understood to follow the r 1 / 4 surface bright- 
ness profile (Simien & de Vaucouleurs 1986). Early spectroscopy (Kormendy & 
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Illingworth 1982) placed large bulges along the oblate rotator line in the rotational 
support V/<Jq vs e diagram (some bulges in barred galaxies were found above the 
oblate-rotator line; Kormendy 1982). The stellar populations of bulges were un- 
derstood as having metallicities higher than solar, after the pioneering work by 
Whitford (1977) and subsequent studies of the MW Baade's Window by Rich and 
collaborators (e.g. Rich 1988). The high metallicities are properties typical of giant 
elliptical galaxies. 

The similarities of bulges and ellipticals strengthened the approach of galaxian 
component analysis whereby a galaxy can be analyzed into largely independent 
structural units, foremost disks and spheroids. Spheroids would encompass bulges 
and ellipticals, sharing similarities in structure, internal dynamics and populations. 
Such similarities would imply similar formation histories and presumably old ages. 

4 Surface brightness profiles 

There is now ample evidence that bulges do not follow the r 1 / 4 law. Kent et al. 
(1991), using Spacelab IR telescope data, found that the vertical profile of the MW 
bulge at 2.5 /im is well fit with an exponential law. Andredakis & Sanders (1994) 
and de Jong (1996) showed that bulge-disk decompositions of external galaxies led 
to more stable fits and smaller residuals when the exponential law was used for 
bulges rather than the r 1//4 law. 

A breakthrough in the characterization of bulge surface brightness profiles came 
with the application of the Sersic (1968) r 1 ' n model, 

J(r) = h ■ exp {-&„ • [(r/re) 1 /" - 1]}, (4.1) 

where r e is the effective radius, I e is the brightness at r e , n is the shape index 
and b n w 1.9992n - 0.3271 (Caon et al. 1993). APB95 showed that If -band bulge 
profiles are well fit with r 1 /™ models, with a continuous distribution of shape 
indices, ranging from n ss 4 — 6 for early types to n « 0.5 for late types (Fig.rjJ. 
Later fitting by G01, Khosroshahi et al. (2000a), Mollenhoff & Heidt (2001), and 
Trujillo et al. (2001) confirmed this dependency. The ability of Sersic's law to fit 
the range of structures of bulges is comparable to its ability to fit elliptical galaxies 
(Caon et al. 1993), suggesting a common range of structures for galaxian ellipsoids 
of a wide range of luminosities. 

Balcells et al. (2003a, hereafter B03) reanalyzed the sample studied by APB95 
using HST/NICMOS F160W images to extend the NIR profiles inward to -0.1 
arcsec. At the increased resolution, Sersic+exponential modeling yields strong 
central residuals: galaxy centers are brighter than the inward extrapolation of the 
Sersic model. B03 show that for most cases (63%) adding a single point source 
(PSF-convolved) provides a good fit (Fig[2). In some cases, a central exponential 
component, instead of a point source, improves the fit. Such central components, 
with colors corresponding to reddened stellar populations (Peletier et al. 1999), 
are largely unresolved by the NICMOS Camera 2 (PSF = 0.19 arcsec FWHM). 
Although the Nuker model (Lauer et al. 1995) can provide good fits to the inner 
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Fig. 1. Points: the distribution of 
bulge Sersic index n against the bulge- 
to-disk luminosity ratio B/D (APB95). 
Arrows: growth vectors for bulges de- 
rived from iV-body simulations of satel- 
lite accretion onto exponential bulges 
(Aguerri et al. 2001). Arrows of increas- 
ing length correspond to M sa t/Mt, u i ge = 
1/6, 1/3, 1, respectively. 



Fig. 2. (a) #-band HST+ground- 
based profile of NGC 5443 (Sb) fitted 
with PSF-convolved Sersic bulge and ex- 
ponential disk components. (b) The 
same profile with a PSF-convolved point 
source (dashed line) added to the fit- 
ting function. Residuals in lower panels. 
From Balcells et al. (2003a). 



regions of the profiles, out to a few arcsec, it cannot model the entire bulge profile. 
We address nuclear components in bulges in §0 

Once nuclear components are included in the fitting function, the best-fit Sersic 
models have considerably lower shape incides n than obtained from ground-based 
profiles. For the B03 galaxies, the n distribution has an upper envelope that rises 
toward more luminous bulges, and a highest value of n = 3.0 is found. 

Bulge structure shows strong links to internal dynamics. Trujillo et al. (2001b) 
define a concentration index, CV c (a) , as the ratio of luminosities inside r e and 
inside ar e (a < 1). For Sersic profiles, C rc (a) is a monotonic function of the 
shape index n. CV e (a) is useful because of its low sensitivity to the depth of the 
images. C Ti! (a) is shown to correlate strongly with the central velocity dispersion 
and with the mass of the Supermassive Black Hole inferred in galaxy centers from 
kinematics measurements (Graham et al. 2001). The correlation with velocity 
dispersion opens an avenue for measuring the Fundamental Plane of spheroids 
using imaging data only (Graham 2002; see also Khosroshahi et al. 2000b). 

4.1 The scaling of bulge and disk parameters 

Studies of well-defined samples of disk galaxies have clarified that the properties 
of bulges are intimately related to those of the disks. While common sense would 
suggest that earlier-type bulges are larger in relation to their disks, it turns out that 
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the ratio of bulge effective radius to disk scale length is independent of Hubble type: 
r e /h w 0.1 (de Jong 1996; Courteau et al. 1996; Graham 2001; MacArthur et al. 
2003). Also, bulge colors are very similar to the inner colors of their disks (Peletier 
& Balcells 1996). These two relations are hard to establish if disks accreted onto 
previously- formed old spheroids; rather, they suggest a coupled formation of bulges 
and disks. 

5 Nuclear structures 

HST imaging at visible and NIR wavelengths has provided an important increase 
in spatial resolution to study galaxy nuclei. At the distances of typical nearby 
samples (cz < 3000 km/s), ground-based images (1" FWHM seeing) resolve 150 
pc, while WFPC2 (PSF 0.1" FWHM) resolves nuclei down to -15 pc. Bulges 
have been targeted in programs by Carollo and collaborators, and by Peletier and 
collaborators, in addition to many programs on ellipticals which also image SO 
galaxies. 

At HST resolution, galaxy nuclei show a wealth of structure that contrasts with 
the conception of bulges as smooth spheroids, including nuclear spirals, dust, star 
formation and central sources often unresolved at HST resolution (Carollo et al. 
1998; Peletier et al. 1999; Ravindranath et al. 2001; Rest et al. 2001; Seigar et al. 
2002). Because the nuclei of ellipticals show similar levels of complexity (Phillips 
et al. 1995; Ravindranath et al. 2001), such structures are often associated with 
the galaxy's bulge. Their vertical sizes are smaller than the disk vertical extent, 
hence we might associate them with the bulge-spheroid or with the inner disk. 

Unresolved sources have been detected in every study of galaxy centers using 
HST imaging (Phillips et al. 1995; Faber et al. 1997; Carollo, Stiavelli & Mack 1998; 
Ravindranath et al. 2001; B03). Figure EH shows point-source absolute magnitudes 
against galaxy absolute magnitude from various sources. Both the frequency of 
detection and the inferred luminosities depend on the method applied to subtract 
the underlying bulge light, so results differ among the studies. The frequency of 
detection ranges from 50% to 85%. Point source luminosities are of order of 10~ 3 
of the bulge luminosities, and, in B03, they scale as Lps — (Lbui) ' 5 - 

6 Kinematics 

Kinematic studies have provided strong evidence that not all bulges follow the 
oblate-spheroid picture derived from V/a vs. e studies in the eighties (Kormendy 
k Illingworth 1982). 

Kormendy (1993) summarizes the evidence for disk- like dynamics in bulges. 
This includes flatness, low velocity dispersion and rapid rotation, as well as spi- 
ral structure. An example of the type of analysis involved is given in Erwin et 
al. (2003) on the bulges of two barred galaxies (NGC 2787, NGC 3945). The 
bulges show flattened light distributions with exponential surface brightness pro- 
files, and high rotation parameters (V/a)* ~ 1.5, which hint that the light previ- 
ously adscribed to their bulges comes instead from luminous disks embedded in 
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Fig. 3. //-band absolute magnitudes 
of nuclear unresolved sources against 
galaxy B-band absolute magnitudes, 
for samples of ellipticals (open circles, 
Ravindranath et al. 2001), SO's (crosses, 
Ravindranath et al. 2001) bulges (filled 
circles, B03). For the latter, bulge 
rather than total magnitudes are plot- 
ted. 
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Fig. 4. Position-velocity diagrams for 
the peanut-shaped galaxy NGC 5746, 
showing the figure-of-eight pattern 
given by gaps between orbit families in 
a bar (Bureau & Freeman 1999). 



the bars. Smaller, non-exponential, rounder structures are found inside the in- 
ner disks. Falcon-Barroso et al. (2003), working on a sample of 18 inclined disk 
galaxies, find that minor-axis velocity dispersion gradients decrease with increas- 
ing bulge ellipticity. Such flat dispersion profiles in high-ellipticity bulges link 
smoothly with dispersions in the disk, indicating a continuity of dynamical tem- 
perature from bulge-dominated to disk-dominated regions of the galaxies. This 
behavior is not restricted to late types; it occurs in galaxies as early as SO and 
with high photometric B/D. 

Boxy or peanut-shaped bulges show a cylindrical velocity field, and very low ve- 
locity dispersions except in a central peak (Shaw et al. 1993), linking box-peanut 
bulges to thickened disks. The connection of these bulges to evolved disks was 
clearly established with the use of high-order kinematic analysis. Kuijken & Mer- 
rificld (1995) and later Bureau & Freeman (1999) applied the Unresolved Gaussian 
Decomposition (UGD) technique to analyze gas and stellar kinematics of edge-on, 
box-peanut bulges, showing the figure-of-eight pattern in the position-velocity di- 
agrams expected from XI and X2 orbits in bars (Fig.^J. That such a feature does 
not appear in bulges with elliptical isophotes strongly links box-peanut bulges to 
bars. 

The above studies reveal that some photometric bulges are really disks with 
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brightness in excess of the inward extrapolation of the outer exponential disk. 
The frequency of kinematically disky bulges can be established by studying well- 
defined samples with the upcoming generation of 3D spectrographs, e.g. SAURON 
(de Zeeuw et al. 2002). 

Bulges rotate around the same axis as the disk does, i.e. they do not show 
minor axis rotation (Falcon-Barroso et al. 2003), with a few notable exceptions. 
As an example, the bulge of NGC 4698 rotates about an axis perpendicular to 
the disk rotation axis (Bertola et al. 1999). These objects indicate that accretion 
or merging has been important for some bulges. The misaligned fraction may be 
related to the frequency of accretion events. However, while monolithic collapse is 
unlikely to lead to disk-bulge kinematic misalignment, mergers do not necessarily 
generate kinematic misalignment. Hence, the fraction of misaligned bulges is a 
lower limit to the fraction of bulges with a rich accretion history. 

6.1 New kinematic spectral analysis techniques 

Falcon-Barroso et al. (2003) have succesfully used single-burst stellar population 
(SSP) models as templates for kinematic spectral analysis. This has become pos- 
sible thanks to the development of new SSP NIR spectral energy distributions 
(SED) at 1.5A resolution (Vazdekis et al. 2003; see also Vazdekis 1999). The 
method consists of measuring the spectrograph's broadening function in calibra- 
tion arc exposures, and convolving the synthetic spectral templates to match that 
broadening. Determining galaxy kinematics can often be done better with SSP 
templates than with stellar templates, probably as a result of reduced template 
mismatch. The method allows a reduction in telescope time requirements, and 
has strong potential for high-redshift galaxy kinematics, where target galaxy and 
stellar templates would be observed in different wavelength ranges. 

7 Fundamental Plane, and bulge population ages 

The distribution of a galaxy sample on the Fundamental Plane (FP) defined by 
old ellipticals can in principle provide three types of diagnostics. The FP is usually 
expressed as 

log(r e ) = cdog(<7o)+/3<M>e+7 (7.1) 
while, assuming Virial equilibrium, one has 

M 

log(r e ) = 2 log(ao) + 0.4 < ft > e + iog(C r C v ) - Iog(— ) + 7' (7.2) 

(Djorgovski, de Carvalho & Han 1988; Bender et al. 1992; Graham & Colless 
1995), where variations of the structural and kinematic parameters C r and C v 
reflect non- homology. First, the thinness of the distribution provides a measure of 
the sample homogeneity. Second, if one believes that all the objects in the sample 
are structurally and dynamically homologous, M/L trends with luminosity may 
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be inferred (leading to M/L ~ L 25 , Faber et al. 1987). This provides population 
age trends with luminosity. Finally, one may investigate whether broken homology 
can reconcile the FP with the VT prediction without recourse to M/L variations 
(Hjorth & Madsen 1995; Capelato et al. 1995). 

Key difficulties in establishing the FP of bulges are the disk contamination to 
the surface brightness and the velocity dispersion, and dust extinction. Falcon- 
Barroso et al. (2002) address these problems by measuring structural parameters 
from fC-band images, and using colors unaffected by dust to convert those to the 
standard B-band. These authors find a small but significant offset of early-to- 
intermediate type bulges with respect to the Coma ellipticals of Jorgensen et al. 
(1996). Their result rests on an exquisite attention to details such as aperture 
corrections for the velocity dispersions. Similar offsets had been found by Bender 
et al. (1993) using B-band imaging data. 

Can the FP offset of bulges w.r.t. Coma ellipticals be due to non-homology 
rather than to population differences? Gonzalez (2003) shows that, for a range of 
ellipsoids produced by collisionless mergers of spirals, C r and C v (ean. l7.2| l do vary 
from object to object, but the product (C' r C' v ) remains constant among merger 
remnants. The structure and internal dynamics conspire to bring non-homologous 
ellipsoids onto the same FP (see also Trujillo, Graham, & Caon 2001), implying 
that merger-driven non-homology may not contribute to offsets or thickening of 
the FP. In terms of kinematic non-homology only, Falcon et al. (2002) find that 
the non inclusion of a rotational support term in the FP can explain the observed 
offset of bulges w.r.t. ellipticals. 

The FP results are generally interpreted as evidence that bulges are, on the 
whole, slightly younger than cluster ellipticals. By how much? When assuming 
homology and ignoring rotation, Falcon et al. (2002) derive bulge ages of around 
10 Gyr if ellipticals are 12.5 Gyr old. This result is consistent with those obtained 
from colors and from line strengths. Optical colors suffer the well-known age- 
metallicity degeneracy. In bulges, dust in the disk only aggravates the problem. 
The degeneracy can however be partially broken by the combined use of optical 
and NIR colors, e.g. U — R vs. R — K or B — I vs. I — H (Bothun & Gregg 
1990; Peletier & Balcells 1996; Peletier et al. 1999). U - R is sensitive to age 
and dust, while R — K is sensitive to dust and metallicity. Fitting single-stellar 
population (SSP) models to colors on dust-free lines of sight, Peletier et al. (1999) 
infer bulges ages above 10 10 Gyr, at most about 2 Gyr younger than the ages of 
Coma ellipticals. Most importantly, no clearly young bulges are found in a sample 
of early- to intermediate- type bulges (Sbc galaxies do show slightly younger ages). 
Line-strength indices yield similar results; in an ongoing program of minor-axis 
spectroscopy of bulges, using H7 to disentangle age from metallicity, Goudfrooij, 
Gorgas, & Jablonka (1999) also find bulge population ages slightly younger than 
cluster ellipticals, and with a similar age span as field ellipticals. 

Dating bulge populations to answer who formed first, the bulge or the disk, is a 
subtle problem. Peletier & Balcells (1996) infer that bulges are older than the inner 
regions of the disks by at most A log £ = 0.11 (Ai = 3 Gyr for tdisk — lOGyr), 
assuming equal metallicities, and by much less if bulge metallicities are higher. 
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One may argue that SSP models applied to local galaxy samples are too simplistic 
to resolve the small age differences between bulges and diks, but the results may 
not be too far off. The ages given above would put the bulk of the star formation 
in bulges and disks at redshifts z > 1.5, a result that is consistent with the onset 
of structurally mature disk galaxies at z ~ 1 (Lilly et al. 1998). A powerful 
avenue to improve these results will be to study epochs closer in time to the epoch 
when the bulge-disk strucures assembled. This demands studying galaxy samples 
approaching z = 1 , something now feasible with recourse to HST optical and NIR 
imaging (e.g. Ellis et al. 2001). 

8 Formation models 

8.1 Accretion 

In a soft version of the merger origin proposed by CDM-based galaxy formation 
models, bulges may form/grow by accreting smaller structural units. Evidence 
for accretion comes from, e.g., orthogonal bulge kinematics (§EJ. Orbital decay 
by dynamical friction will indeed feed the galaxy's central bulge. The question is 
whether the numbers add up: i.e. did satellites exist with the right numbers and 
the right masses for dynamical friction orbital decay to build up a bulge mass in 
a Hubble time?. Fu et al. (2003) investigate orbital decay of Super Star Clusters 
(SSC) to explain the formation of small, late-type bulges with cuspy cores. At the 
other end of the B /D spectrum, accretion formation models for large bulges were 
studied by Pfenniger (1993). He modeled the accretion of several satellites onto 
a disk-halo galaxy. The satellites end up at the core, while the disk thickens to 
a Sombrero-like shape. His model, with rigid satellites which do not lose mass, is 
probably simplistic; nevertheless, it highlights that bulges grown by accretion are 
not composed of stars from the accreted galaxies only, and that other mechanisms 
besides bar buckling instabilities may pump disk material up into a bulge-like 
distribution. 

Aguerri, Balcells, & Peletier (2001) investigate the accretion hypothesis by 
checking whether bulges grown by accreting satellites follow the trend of Sersic 
index n with B/D observed by APB95. They study whether the accretion of 
satellites onto model galaxies harboring exponential bulges provides not only ad- 
ditional mass to the bulge but also a structural change toward higher-n Sersic 
profiles. The results of the TV-body models are shown in Figure ^ with arrows 
tracing growth vectors for bulges in the n-B/D plane, together with the points 
from APB95. It is promising that the trends are the same, although growth of 
n appears steeper in the models. As it turns out, exponential profiles are very 
fragile to merging. Hence, it is unlikely that bulges with exponential profiles grew 
significantly through collisionless mergers. The models have several limitations. 
The disks thicken during satellite accretion, and rebuilding of a thin disk needs to 
be postulated (the thick disk left over after the accretion may be related to the 
thick disks of early-type disk galaxies) . It is not clear how bulge growth via merg- 
ers plus rebuilding of the thin disk can explain the observed scaling of bulge and 
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disk parameters (S I4.1|I . Accretion of satellites could explain the class of distorted, 
thick boxy bulges studied by Liitticke et al. (2000). 

8.2 Coupled bulge-disk formation models 

Van den Bosch (1998, 2000) included bulges in the Fall & Efstathiou (1980) theory 
of disk formation. Gas cooling inside a virialized halo is highly unstable and 
forms a central bulge component rather than a disk. When the bulge has grown 
sufficiently, its central potential allows subsequent gas to settle into an ordered 
disk. The stability condition determines a higher bulge-to-disk ratio at increasing 
formation redshift. This model broadly accounts for the ages of bulges and disks, 
and allows for the formation of a central bulging component without recourse to 
merging; the bulge stars did not dominantly form in disks as they do in the merger 
hypothesis. 

9 Conclusions 

The ability of the Sersic model to describe the shapes of the surface brightness 
profiles of both bulges and ellipticals suggests that bulges and ellipticals belong 
to a single family of ellipsoidal galaxian components. The proximity of the FP 
of both bulges and ellipticals also emphasizes their similar structure. Merger 
models account for some of these properties. On the other hand, both the detailed 
structure and kinematics of some bulges of all Hubble types display properties 
more typical of thickened disks than of oblate spheroids. Kinematic studies of 
large, well-defined samples will be useful for determining the prevalence of disky 
"bulges" across the Hubble sequence and the bulge luminosity sequence. Once we 
know how to distinguish disky from classical spheroidal bulges, a new term might 
be desirable to refer to disky bulges. 

Structural parameters and colors of bulges and disks are coupled, arguing 
against a clean separation of the epochs and formation mechanisms for bulges 
and disks. Population diagnostics show that populations in early-type bulges are 
old, with ages above 10 10 Gyr, approaching those of cluster ellipticals. These two 
results place the assembly of disk-bulge galaxies at redshifts z > 1.5. 
Acknowledgements: I would like to thank R. Peletier, A. Graham and P. Erwin 
for discussions and for comments on the manuscript. 
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